Do shape invariant solitons in highly nonlocal nematic liquid crystals really exist? 
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We question physical existence of shape invariant solitons in three dimensional nematic liquid crys- 
tals. Using modified Petviashvili's method for finding eigenvalues and eigenfunctions, we determine 
shape invariant solitons in a realistic physical model that includes the highly nonlocal nature of 
the liquid crystal system. We check the stability of such solutions by propagating them for long 
distances. We establish that any noise added to the medium or to the fundamental solitons induces 
them to breathe, rendering them practically unobservable. 
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The fundamental spatial optical soliton is a beam that 
propagates in a nonlinear (NL) medium without changing 
its transverse profile [1]. Such shape-invariant solutions 
are easily identified in (l+l)-dimensional [(1+1)D] NL 
systems because the inverse scattering theory [2] guar- 
anties their existence. The situation is less clear in the 
multidimensional and multicomponent systems. No cred- 
ible inverse scattering theory is formulated in more than 
one dimension and even when the localized solutions are 
found, no credible procedure for guaranteeing their sta- 
bility is established (however, there are many credible 
linear stability analyses or stability criteria or numerical 
procedures, but not for rigorously proven stability.) In 
fact, wave instability and the collapse of solutions are 
overriding concerns in multidimensional NL systems 
Additional compounding difficulties arise in the multi- 
component vector models or in the scalar nonlocal mod- 
els in which the medium response is driven by the optical 
field itself. Such are the models describing the genera- 
tion of solitary waves - nematicons - in nematic liquid 
crystals (NLCs). 

Nonlocality is an important characteristic of many NL 
media. A highly nonlocal situation arises in a nonlocal 
nonlinear (NN) medium when the characteristic size of 
the response is much wider than the size of the excitation 
itself 0, H[ . In NLCs both experiments 0, [?J and theoret- 
ical calculations d, Q demonstrated that the nonlinearity 
is highly nonlocal. 

For more complex nonlinear it ies, numerical techniques 
are necessary to determine the soliton solutions. Soliton 
profile calculations in NN media have been presented in a 
number of papers, see for example flol-[l2j. The existence 
and stability of 2D solitons in media with NN response 
was discussed in [13]; even a high degree of nonlocality 
did not guarantee the existence of stable high-order soli- 
ton structures [l4|. Orientational nonlinearity in NLCs 
is highly nonlocal but the NL response is not perfectly 
quadratic, implying that if one launches a Gaussian beam 
into the cell it is only possible to observe breathing soli- 
tons 0,113. 



In some publications soliton profiles were calculated 
using semi-analytical models 0, fl6l-[l8j. For the more 
general vectorial model, in which the order parameter in 
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FIG. 1: Liquid-crystal cell model adopted. 



NLC is not constant, steady elliptical soliton profiles are 
found numerically in [19j. To determine such profiles, the 
authors demonstrated that it is necessary to include all 
three components of the optical electric field. 

However, what is puzzling is that even though ev- 
erybody agrees that shape-preserving solitons do exist 
in highly nonlocal NLCs, practically nobody cared to 
present them explicitly. Experimental accounts profusely 
mention steady nematicons, but careful inspection of all 
published figures reveals self- focusing oscillations. True, 
experimental results may be of not much help in this re- 
gard, because all experimental setups feature a few mm 
long cells, which cannot capture slow (if any) conver- 
gence to a steady profile. In this paper we find a family 
of fundamental solitons for the same model and the same 
parameters; we check their stability in propagation and 
demonstrate that any small change in the input shape, 
as well as in the medium, leads to the soliton breathing. 
Consequently, we question the real physical observability 
of such shape- invariant solitons. 

To find an exact fundamental soliton solution in a NLC 
model that is not vectorial, we use an iterative numerical 
eigenvalue technique. We consider widely accepted scalar 
model of beam propagation that is well suited for uniaxial 
NLCs and low-intensity optical fields, which correspond 
to most situations of practical interest. Also, we discuss 
the influence of boundary conditions (BCs) on the shape 
and power of solutions, and analyze soliton and Gaussian 
propagation using two different propagation methods. 

We adopt the well-known NN 3D scalar model in 
NLCs, which provides good agreement with experimen- 
tal data Q • The optical beam polarized along the x axis 
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FIG. 2: Fundamental soliton solutions for \i — 5L^ X . Inten- 
sity profiles and reorientation angle distributions are shown 
for two different BCs. Other parameters are the same. 



propagates in the z direction, while the NLC molecules 
can rotate in the x-z plane. The liquid-crystal cell of 
interest is sketched in Fig. [TJ The total orientation 
of molecules with respect to the z axis is denoted as 
0(x,y,z), whereas the orientation induced by the static 
electric field only is denoted by 9q (the pre-tilt angle). 
The bias field points in the x direction and is uniform in 
the z direction; hence the pre-tilt angle is uniform along 
the z axis as well. The quantity = — 0o corresponds 
to the optically induced molecular reorientation. 

The system of equations of interest consists of the 
scaled NL Schrodinger-like equation for the propagation 
of the optical field A, and the diffusion equation for the 
molecular orientation angle @, H, [H| : 



dA 

2i— + A^A + a[sin 2 #-sin 2 #oM = 0, 
oz 

2A^#+[/3 + a|A| 2 ]sin(2#) = 0, 



(1) 
(2) 



where the coefficients a and f3 are proportional to 
the optical and static permittivity anisotropics of the 
NLC molecules, respectively. Hard boundary con- 
ditions (BCs) on the molecular orientation at the 
NLC cell faces in the x direction are assumed: 
6{x = -D/2,y) = 6(x = D/2,y) = const. [20]. In our 
calculations we use data corresponding to typical exper- 
imental conditions [H, [TH, Hq| • 

The solitary eigen-solutions are determined from the 
system of Eqs. (|l|2p using the modified Petviashvili's 
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FIG. 3: (a) Fundamental soliton intensity profile obtained by 
the eigenvalue method (black dots), fitted with a Gaussian. 
Parameters: P = 10.6 mW, fi = 3.84L^; zero BCs. (b) Soli- 
ton and Gaussian propagation using two different propaga- 
tion methods, FFT+SOR and FFT+FFT. The soliton power 
P 10.6 mW; Gaussian power P=10.6 mW for FFT+SOR, 
P=10.1 mW for FFT+FFT. 



iteration method [2ll-[23|. Equation (pQ) suggests the 
existence of a fundamental soliton of the form A — 
a(x,y)e lflz , where /i is the propagation constant. The 
real-valued function a(x, y) satisfies the equation: — Aa + 
(2/i + P)a = Q, where P = asin 2 ((9 ) and Q = a sin 2 (0)a. 
After Fourier transforming the equation for a, we get: 



■2/i 



(Q-Pa), 



(3) 



where overbar denotes Fourier transform. Straightfor- 
ward iteration of Eq. ((3]) does not converge in general, 
so the stabilizing factors had to be introduced [22|, [23| . 
In each iteration step of Eq. (|3j), Eq. (|2j) is treated using 
a successive overrelaxation (SOR) method, until conver- 
gence is achieved. 

Stable soliton solutions are presented in Fig. [2] for two 
different BCs. The shape and the power of fundamental 
shape-invariant solutions naturally depend on the BCs 
applied. Zero BCs (0 = on all boundaries) correspond 
to the Dirichlet BCs. Periodic BCs correspond to the 
mixed BCs - Dirichlet along the y axis and Neumann 
along the x axis. The solution with the periodic BCs 
is more appropriate to the geometry of the problem; fur- 
thermore, it is more acceptable on physical grounds. The 
fundamental soliton so obtained requires less beam power 
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FIG. 4: Propagation of the fundamental soliton in a noisy 
medium. An amount of 0.5% randomly changing noise is 
added to 0q. Red sinusoidal fit is to guide eye. 

for the same value of the propagation constant and iden- 
tical other parameters. In identical conditions (but for 
BCs) the solution requiring less power should be favored. 
The often used 0q = 7r/4 approximation leads to the so- 
lution with zero BCs, and consequently such a soliton is 
less appropriate. 

Spatial solitons in highly nonlocal media with 
quadratic response possess Gaussian profiles [H, H| . How- 
ever, the fundamental soliton profile is not Gaussian. The 
soliton intensity profile compared to a Gaussian is shown 
in Fig. [3ja); the difference is confined to the tails. To 
check the stability of fundamental solitons, we propagate 
them numerically; the results are presented in Fig. [3fb). 
Also included in Fig. [3fb) is a case presenting propaga- 
tion of a Gaussian with similar parameters, but obtained 
using two different numerical methods. In both methods 
a split-step beam propagation procedure based on the 
fast Fourier transform (FFT) is used for the propaga- 
tion of the optical field. In the first method the diffusion 
equation for the optically induced molecular reorienta- 
tion is treated using the SOR method. In the second 
method the diffusion equation is treated using the split- 
step procedure again. One can see that the methods 
provide similar results; however the first method is more 
accurate. 

The problem with the FFT procedure is that it treats 
an array of transversely periodic cells. Since the molec- 
ular reorientation is wide, it tends to slightly spill over 
into the adjacent cells, i.e. back onto itself, adding to the 
optical field. This is not an overriding problem in the 
propagation of a Gaussian, as it only leads to a slightly 
amplified oscillation of the breathing solution. However, 
it makes huge difference in the propagation of the fun- 



damental soliton - it makes it impossible for the field 
to keep the shape-invariant input profile and therefore 
should be discarded. Even the SOR solution slightly os- 
cillates at lower accuracy; this, however, becomes imper- 
ceptible as the accuracy is improved. In Fig. [3fb) we 
show a case where the oscillation of the amplitude is still 
perceptible. This brings us to an important point. 

When one considers the propagation of a Gaussian 
beam using the two propagation methods, the results are 
close. The propagation of Gaussians invariably leads to 
breathing beams, regardless of the method of integra- 
tion. By the same token, when the fundamental soliton 
is propagated through the medium in which a small ran- 
dom noise added to the underlying molecular orientation 
#o, a breathing solution is also obtained. The same phe- 
nomenon happens as well when a small intensity noise is 
added to the fundamental profile, but 0q kept unchanged. 
This phenomenon is confirmed in our computations (Fig. 
2J and is not difficult to understand. In a highly NN 
medium any additional energy from noise, no matter how 
small, cannot be radiated away and the solution has no 
way to relax to the fundamental soliton. Therefore, it 
keeps oscillating about the fundamental soliton, forming 
a stable breathing soliton. Since noise is unavoidable in 
any realistic set-up, be it experimental or numerical, this 
fact opens the question of the physical observability of 
shape-invariant fundamental solitons in highly NN me- 
dia. 

In conclusion, we have presented calculations of shape- 
invariant fundamental solitons in a highly nonlocal 3D 
scalar NCLs, for a realistic physical model. Using modi- 
fied Petviashvili's iterative scheme we numerically deter- 
mined the fundamental spatial soliton profiles and found 
a family of solutions, depending on BCs. We depicted 
stable propagation of such solitons. We then demon- 
strated that upon propagation, any amount of noise 
transforms the fundamental soliton into a breather, mak- 
ing shape-invariant nematicons practically unobservable. 
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